Fractional topological superconductors with fractionalized Majorana fermions 
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In this paper, we introduce a two-dimensional fractional topological superconductor (FTSC) as a 
strongly correlated topological state which can be achieved by inducing superconductivity into an 
Abelian fractional quantum Hall state, through the proximity effect. When the proximity coupling 
is weak, the FTSC has the same topological order as its parent state and is thus Abelian. However, 
upon increasing the proximity coupling, the bulk gap of such an Abelian FTSC closes and reopens 
resulting in a new topological order: a non-Abelian FTSC. Using several arguments we will conjec- 
ture that the conformal field theory (CFT) that describes the edge state of the non-Abelian FTSC 
is U{l)/Zi orbifold theory and use this to write down the ground-state wave function. Further, 
we predict FTSC based on the Laughlin state at v = 1/m filling to host fractionalized Majorana 
zero modes bound to superconducting vortices. These zero modes are non-Abelian quasiparticles 
which is evident in their quantum dimension of d m = y2m. Using the multi-quasi-particle wave 
function based on the edge CFT, we derive the projective braid matrix for the zero modes. Finally, 
the connection between the non-Abelian FTSCs and the Z^m rotor model with a similar topological 
order is illustrated. 



INTRODUCTION 



There is growing interest in strongly correlated topo- 
logical states of matter and their potential application 
in topological quantum computation. [TI-fTT] Interaction 
can lead to novel topological states exhibiting fractional 
excitations, anyon statistics of quasiparticles, and non- 
trivial ground-state degeneracy (GSD).[T2J [T3] Interest- 
ingly, these properties of the topological states arc ro- 
bust and insensitive to disorder or interaction as long as 
the bulk gap does not close. [T] An appealing topological 
state of matter is the non-Abelian phase that supports 
anyon excitations with quantum dimensions more than 
one that obey non-Abelian braid statistics. |14ftl6] This 
exciting possibility along with the robustness of the topo- 
logical states that immunizes anyons from quantum de- 
coherence make non-Abelian anyons a promising candi- 
date for topological quantum computation [1 . The first 
non-Abelian topological phase was proposed by Read and 
Moore for the fractional quantum Hall (FQH) state at 
v = 1/2 filling fraction. Their so-called Pfaffian wave 
vfunction is a paired state of composite fermions whose 
excitations obey Ising-type non-Abelian statistics. [T3] 
It was shown later that the projective Ising-type non- 
Abelian statistics can also be realized in the p + ip su- 
perconductors in the weak pairing phase. [T7J [TB] In- 
deed, the Pfaffian FQH state and a p + ip superconduc- 
tor have the same topological order, whose excitations 
are the celebrated Majorana fermions with quantum di- 
mensions equal to 

Recently, it has been shown that the p + ip supercon- 
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FIG. 1: A schematic phase diagram of the FTSCs. gi is a 
tuning parameter that drives the phase transition between an 
FQH at v = 1/m filling fraction and a trivial insulating state, 
e.g., a Wigner crystal. 32 denotes the tunneling amplitude 
between the superconducting substrate and the FQH sample. 
The edge theory of the Abelian FTSC side is described by 
a free boson compactified on a circle with radius R = \/m, 
i.e., U(l)m CFT, while the edge state of the quantum critical 
region i.e. the non-Abelian side, is given by the U(V) m /Z% 
orbifold CFT. The trivial superconductor does not support 
gapless edge states. 



ductor can be achieved by inducing superconductivity in 
the bulk of an integer quantum anomalous Hall (QAH) 
state via proximity to a superconductor. 19, 2U] In a dif- 
ferent study, Bombin showed that Majorana fermions can 
emerge in an Abelian topological state with topologi- 
cal defects. [21] These two mechanisms can be related 
through viewing superconducting vortices as twist oper- 
ators. These observations motivated us to investigate the 
fate of an Abelian FQH state or a fractional Chern insu- 
lator at v = 1/m filling fraction in a similar setup. It will 
be shown that the weak tunneling of Cooper pairs does 
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not affect the topologicaf nature of the FQH state. We re- 
fer to the resulting phase as the Abelian fractional topo- 
logical superconductor (FTSC). However, we discuss that 
the system can undergo a topological phase transition by 
increasing the tunneling strength. The new topological 
state will be a non-Abelian FTSC, whose excitations are 
fractionalized Majorana fermion zero modes, «o- These 
non-Abelian excitations are bound to the superconduct- 
ing vortices, follow ccq = 1 relation, and have quantum 
dimension d m = \2m. 

The paper is organized as follows. In Sec. II, we define 
the notion of a FTSC and try to build a model for it. 
To make the situation more comprehensible, an exactly 
solvable model will be reviewed in which a topological 
phase transition happens between an integer QAH and 
a p + ip superconductor. In Sec. Ill, we identify the 
CFT that describes the edge theory of the FTSCs. The 
correspondence between the bulk wave function and the 
conformal blocks of the edge theory is explained, and the 
notion of the twist field is illustrated with an emphasis on 
their crucial role in the topological order. In Sec. IV, we 
utilize the CFT that describes the non-Abelian FTSCs 
to study its topological order. The existence of the zero 
modes bound to vortices of the superconducting state is 
proven. Furthermore, the quantum dimension and the 
braid statistics of the zero modes is derived. Sec. V is 
devoted to the computation of the bulk wave function of 
the non-Abelian FTSCs. In Sec. VI, we study the gauge 
theory of the non-Abelian FTSCs. It turns out to be a 
Zim gauge theory which allows us to make a connection 
between our prediction about the quantum dimension of 
the fractionalized Majorana zero modes and the results 
of the recently studied rotor model with a similar gauge 
structure. 



II. WHAT ARE FRACTIONAL TOPOLOGICAL 
SUPERCONDUCTORS? 

In this section, we try to present a definition of a FTSC 
and what exactly we mean by this term. Roughly speak- 
ing, an FTSC is a superconducting state of matter whose 
parent state is an FQH state. Two topologically dis- 
tinct classes of FTSCs can be imagined, Abelian and 
non-Abelian FTSCs. Abelian FTSCs are the counter- 
part of the Abelian FQH states and refer to the case 
where all the excitations of the superconducting ground 
state are Abelian anyons. In other words, exchanging 
any of two quasiparticles can only change the phase 
of the ground-state wave function, and two successive 
exchanges are commuting. Accordingly, all quasiparti- 
cles have quantum dimensions equal to one. An impor- 
tant signature of the Abelian FTSCs is the uniqueness 
of their ground-state wave functions in the presence of 
vortices. On the other hand, non-Abelian FTSCs have 
non-Abelian excitations with non-Abelian statistics. The 
quantum dimension of a non-Abelian quasiparticle (non- 
Abelion) is necessarily more than one leading to degen- 



erate ground states. Exchanging two non-Abelions not 
only can change the phase of the ground state, it may 
also act on the space of ground states by a unitary trans- 
formation. In other words, the ground state wc started 
with can superpose with other degenerate ground states 
after exchanging two non-Abelions. 

Based on what we said in the above paragraph, an 
Abelian FTSC has the same topological order as an 
Abelian FQH state. In other words, if we break all the 
symmetries of both phases, they are smoothly connected; 
i.e. one of them can evolve to the other one adiabatically 
without passing any phase transition. An FTSC has all 
the symmetries of its parent FQH state , except the par- 
ticle number conservation [U(l) symmetry]. The U(l) 
symmetry can be easily broken by inducing supercon- 
ductivity in the FQH state either intrinsically or extrin- 
sically. In the latter scenario, the FQH system can be put 
above a superconductor, so that Cooper pairs can tunnel 
into the FQH state through proximity effect. If we keep 
the many-body gap open in the whole procedure, the re- 
sulting state will be an Abelian FTSC. Consequently, we 
expect the same kind of chiral edge theory for the system 
as that of the parent FQH, except that Abelian anyons do 
not carry a definite electric charge anymore. A more in- 
teresting situation can happen when the many-body gap 
closes due to the proximity and the system undergoes a 
topological phase transition accordingly. The important 
question that we are going to address in this paper is the 
fate of the Abelian state. In the absence of the super- 
conductivity, the fate of an FQH is more or less known. 
For instance, upon decreasing the filling fraction, e.g., 
through tuning the chemical potential, we can end up 
the Wigner crystal phase at small filling fractions. Now, 
let us consider the quantum phase transition between the 
FQH state and the trivial insulating phase which can be a 
Wigner crystal, for example. A weak proximity coupling 
turns the two sides of the phase transition an Abelian 
FTSC and a trivial superconductor with no edge state, 
respectively. However, the strong tunneling of Cooper 
pairs can open a quantum critical region between these 
two regions. As is shown later, we conjecture through dif- 
ferent arguments and analogies that the intervening state 
is an FTSC with non-Abelian excitation. It is shown that 
when inducing superconductivity in the Laughlin state 
at v = 1/m filling, the vortices host non-Abelian frac- 
tionalized Majorana zero modes with quantum dimension 
d m = y/2m. 



A. Superconductivity proximity effect on integer 
quantum Hall states 

In this section, we discuss the superconductivity prox- 
imity effect in a chiral model of the noninteracting elec- 
trons whose phase diagram resembles the physics we are 
pursuing. This simple model has been proposed by Qi 
et al. [2D] A similar model has also been studied by Fu 
and Kane.[19 Before writing the Hamiltonian, we would 
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like to point out some comments that makes the results 
more comprehensible. Imagine an integer QAH system, 
e.g., a band structure with Chern number C = 1. The 
edge theory of this system can be described by chiral 
electrons (or free bosons) propagating along the edges 
of the system and localized in the normal direction. By 
tuning the mass term (in this model the magnetization 
of electrons), we can undergo a phase transition into a 
trivial state with C — 0. At the critical point, the system 
becomes gapless and therefore, an electron at the upper 
edge, which is left mover, can hybridize with (or tunnel 
into) the right mover lower edge state. This can happen 
since the bulk state is gapless at the critical point and 
edge electrons can hybridize with the bulk electrons to 
travel from one side to the other. This coupling between 
the upper and lower edges can destabilize the edge state 
and gap it out. When gap reopens by increasing chemi- 
cal potential, the system becomes a trivial phase with no 
gapless edge state. 

Now, let us induce superconductivity in the whole sys- 
tem, i.e. in the bulk as well as the edge. In the weak 
coupling regime where only a few Cooper pairs can tun- 
nel into the system and cannot affect the band-gap, the 
situation is clear. It turns the QAH state into an Abelian 
topological superconductor whose edge state is still de- 
scribed by chiral complex fermions (Bogoliubov quasipar- 
ticles) . Complex fermions, e.g. electrons have two degree 
of freedom, the real and the imaginary part of fermion 
operator. Each of these degrees of freedom can be repre- 
sented by a Majorana operator, which is a real fermion. 
By computing the Chern number of the resulting Abelian 
superconductor, we would obtain C = 2 referring to the 
fact that the edge state can be described by two Majo- 
rana fermions or equivalently a single complex fermion, 
in this case Bogoliubov quasiparticles. Moreover, induc- 
ing superconductivity into the trivial phase via proximity 
effect simply turns it to a trivial superconductor. Now let 
us discuss what happens in the strong coupling between 
the superconductor and the QAH system. By strong we 
mean that the induced Cooper pairs can close and reopen 
up the band-gap in a new fashion. The best way to think 
about this problem is to consider the edge theory when 
the bulk gap closes. Like before, left-mover edge elec- 
trons at the upper boundary can tunnel into the lower 
edge state with right-mover electrons through the gapless 
bulk and hybridization with bulk electrons that obviously 
causes instability. However, note that we have a differ- 
ent instability channel, the Cooper pairing, as well. The 
upper left mover edge electrons can pair with the right 
mover edge electron, since the gapless bulk is a super- 
conducting state and can mediate Cooper pairs. It is 
this Cooper pairing channel that provides novel possibil- 
ities. Now let us look more carefully. Approaching from 
the Abelian topological superconducting region, the edge 
theory is described by a single complex fermion or equiva- 
lently two real (Majorana) fermions. Through the pairing 
between the upper and lower edges, we can gap out one of 
the real (Majorana) fermions and keep the other gapless. 



Therefore, we end up with a topological superconductor 
whose edge theory is given by a single Majorana fermion 
mode. Since Majorana fermions obey non- Abelian statis- 
tics and have quantum dimension d = v2, the interven- 
ing state is a non-Abelian integer superconductor. Now, 
let us present the model Hamiltonian which is given by 
the following equation: 

HQAH = ^2 = ^phQAH (p) ip P , 
v 

hQAH(p) -{A(p x + i Py ) -m(p) J' (1) 

where m (p) = mo + Bp 2 and tp p — (c p ,^,c p ^) T . This 
type of coupling between spin-up and spin-down electrons 
can be achieved, for example, in the Rashba-type spin- 
orbit coupled system, e.g., at the surface of the three- 
dimensional topological insulators, mo represents the 
average magnetization of electrons that can be realized 
either through the Zeeman-type interaction with an ex- 
ternal magnetic field or the exchange interaction with a 
ferromagnet on the top. Assuming all the negative en- 
ergy states are occupied and |m | > 0, the Hall conduc- 
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tivity of the above Hamiltonian is given by o~h = C^, 
where C is the Chern number of the lower energy band. 
It can be shown that C = +1 when mo < and vanishes 
for mo > 0. Now we turn the above QAH system to a 
TSC. An effective pairing interaction between quasipar- 
ticles can be induced in the proximity of an s-wave su- 
perconductor. The corresponding Bogoliubovde Genncs 
(BdG) Hamiltonian is 

TT _ 1 V" V&t ( H Q AH W ' V iA(jV \ M, 



where \t p — ( c Pj f, Cp^, c_ p ^, c_ p j) and the chemical 

potential /i has been introduced so that we can dope 
the unpaired QAH state. Intuitively, we expect that the 
ground state becomes a topological superconductor for 
|A| > | mo | and /i — 0. The full phase diagram of the 
above Hamiltonian has been derived in Ref. [20]. There 
are several interesting facts about them, which we men- 
tion below. 



/i = case 

(A) For A = 0, there is a phase transition at mo = 
0. For m < 0, the edge state hosts a chiral Luttinger 
liquid. This edge state which has central charge c — 1, 
can be written in terms of two chiral Majorana modes 
each having c = §. However, for mo > there is no 
edge state. Therefore, both Majorana edge states gap 
out through the transition at fi = A = 0. 
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(B) For nonzero values of A, when A < —to, we still 
have two gapless Majorana modes at the edge. The wave 
function of these Majorana modes is extended along the 
edge and localized in the perpendicular direction. As 
A approaches —to from below, the localization length of 
one Majorana mode normal to the edge becomes larger 
and diverges at A = — m . However, the other Majorana 
mode is still localized normal to the edge, and extended 
along it. Hence, at A = —too, only one of the Majo- 
rana modes delocalizes and merges with the Bulk states. 
When — to < A < to , that Majorana mode becomes 
gapped and only the other zero mode survives. As we ap- 
proach A = too, the localization length for the remaining 
Majorana mode normal to the edge increases, and finally 
diverges at A — to . For to > A, the only remaining 
Majorana modes gaps out as well. Therefore, a topologi- 
cal superconductor emerges at the critical quantum point 
between a QAH state and a trivial insulator. 

(C) The above transition can be described in the fol- 
lowing way. In the absence of pairing potential (A = 0), 
both Majorana modes have the same correlation length 
perpendicular to the edge (we estimate the width of a 
Majorana mode by this length scale) which is £ oc t^t- 
As a result, both will be gapped at the same critical 
point. However, turning on a nonzero A, the correla- 
tion lengths of two Majorana modes differ and we have 
£1 oc 1/ |m + A| and £2 oc 1/ |to — A|. At A = — to , 
the first one has an infinite correlation length and delocal- 
izes, while the other one has a finite correlation length. 
The first Majorana zero mode merges with bulk states 
and disappears beyond that point. At A = +too, the 
remaining Majorana mode acquires infinite correlation 
length and gaps out accordingly. 

(D) A similar phase diagram is expected for the prox- 
imity of an integer quantum Hall (QH) systems (filled 
Landau levels) to an s-wave superconductor. The only 
topological index that labels the integer QAH or QH 
states is their Chern number; therefore, they are topolog- 
ically indistinguishable. The same edge theory describes 
both, so we expect the same topological phase transitions 
in both cases. 



> case 

In Ref. [50], a phase diagram has been presented for 
this case. However, we argue that their result is valid 
only in the weakly correlated systems where either band 
structure is far from being flat or interaction is weak. 
Going back to the result of Qi et al, at A = 0, nonzero 
chemical potential yields a metallic state. However, turn- 
ing on the A term, chemical potential changes only the 
phase transition point given by the following expression: 

i 3 + /l J = m 2 . (3) 




FQH + SC <> 



I 1 > I 

2,R 

FIG. 2: The quantum critical point (QCP) between an 
Abelian and a non-Abelian FTSCs whose parents are 1/m 
Laughlin state. At the QCP, the bulk gap vanishes and 
the localization length of the edge state diverges. As a re- 
sult, the top (1,L) and bottom (2, 7?) sides of the edge state 
that are counterpropagating can interact in the Cooper chan- 
nel represented by the dashed line. Approaching the QCP 
from the Abelian side, the edge quasiparticles are complex 
anyons with two degrees of freedom: Vi = e 1 ^^™, where 
I £ Zra- This pairing between (1,L) and (2, 72) can gap out 

the (Vi -Vty /2i part of the anyon field. Accordingly, only 

the remaining part, ^V; + V"/j /2 = cos (lip/^/m), which is in- 
variant under the Z2 action, ip — > —tp, remains gapless. This 
argument suggests that the edge theory of the resulting state 
after gap reopens is achieved by taking cos (^/rrup) as the elec- 
tron operator instead of e 1 ^™'*'. The CFT that contains this 
electron operator, and with central charge c = 1, is known to 
be a U(l) m /Z2 orbifold theory. [22l 123] Interestingly, such an 
orbifold CFT contains non-Abelian primary fields with quan- 
tum dimensions \Jlm. 



B. Generalization to the fractional case 

Let us look more carefully at the doped QAH state 
in the proximity of an s-wave superconductor, discussed 
above. Assuming a nearly flat lowest band for the Hamil- 
tonian describing the QAH state Hqah and doping this 
system so that the filling fraction becomes fractional, e.g., 
v = — , we can add an interaction whose coupling con- 
stant is significantly larger than the bandwidth and sig- 
nificantly smaller than the energy gap between lower and 
upper bands. This is the definition of a nearly flat Chern 
insulator. Although, doping leads to a metallic behav- 
ior in the integer Chern insulators, in the present case 
the many body wave function becomes nontrivial and ex- 
hibits topological order. The excitations carry fractional 
charge e* = v at v = ^- filling fraction, and they have 
fractional statistics with 9 = — . Finally, the GSD on the 
torus is to. When the system has an open bound- 

ary, the edge state is described by a chiral boson mode 
ip (x,i), from which we can construct the electron oper- 
ator $ (x,t) = e l V^v(x,t)^ Similarly^ th e anyon operator 
with electric charge ±e^j can be represented in terms 
of the free chiral boson as V±i — e^^"^'^. These 
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complex anyon operators have two degrees of freedom, 
cos -J^<p (x, t) and sin -^<-P (x, t), which are eigenstates 
of the following Z 2 symmetry operation: 1Z : <f> — > —(f). 
Interestingly, these primary fields do not have a definite 
electric charge and are suitable when the electromagnetic 
gauge symmetry is broken, i.e. in a superconductor. The 
key assumption which is used extensively in the rest of 
this paper is that above a critical pairing strength, only 
one of these Zi eigenstates survives and the other be- 
comes gapped. Therefore, the electron operator at the 
edge is V e = = cos (y/mip), which is self-conjugate, 

i.e., V e x V e ~ 1. The CFT that contains this electron 
operator at c = 1 has been studied extensively in the 
CFT literature and is called U(l) m /Z2 orbifold theory. 
[1 [52 [53 I5UH55] This orbifold is known to have non- 
Abelian primary fields, suggesting the name non-Abelian 
FTSC for the bulk state. It is worth mentioning that the 
non-Abelian FTSC can also be obtained by inducing su- 
perconductivity in the FQH states, and the procedure 
explained above is not specific to fractional QAH sys- 
tems. 



III. CFT APPROACH TO THE FTSCs 



the fusion algebra of that CFT, and the bulk wave- 
function. We also study the non-Abelian statistics of 
the fractionalized Majorana zero modes and compute the 
corresponding braid matrix. 

B. Edge theory of the Abelian FTSCs 

As we mentioned before, an Abelian FTSC is smoothly 
connected to its parent FQH state at the same filling 
fraction. Thus, they belong to the same topological class 
and have identical edge theories, except that the number 
conservation is broken at the edge of an FTSC. So it is 
enough to consider the edge CFT of an FQH state. The 
edge of a Laughlin state at v = 1/m is described by a 
chiral gapless mode given by the following action (T2"]: 

m f 

S = -— dxdt (d t + vpd x ) ip (x, t) d x <p (x, t) . (5) 

The above edge theory can be written in terms of a U(l) 
CFT on the 2D space-time. To this end, we use the 
following notations: 



A. Bulk wave function and correlation functions of 
the edge CFT 

For the chiral topological states, it is conjectured that 
the holomorphic part of the bulk wave function can be 
described as the correlation function of the fermion op- 
erators in the chiral edge theory. [TH [TBI [33] First of all, 
we have to identify the CFT that describes the chiral 
edge state. Second, we should find the electron opera- 
tor V e (z) that has a half integer conformal dimension. 
This is to makes sure that two electron operators anti- 
commute. Then the wave function of the ground state 
of N electrons with electrons coordinates {z/,}, where 
zr = Xk + iyki is simply given in terms of the correla- 
tion function involving N electrons operators at {zk} so 

'F ({ z k}) = ^IlfcLi ^ e ( z k)\ Third, we can construct the 
excited wave function by inserting other primary fields in 
the correlation wave function. For example, in the 1/m 
Laughlin states, the electron operator is V e = e l ^™- v and 
V a = e v 7 ™ is the primary field associated with the frac- 
tional excitation whose charge is q — e/m. Hence, the 
holomorphic part of the wave function for N electrons 
and n anyons would be of the following form: 

In N \ 

* dm}; W}) = (j[v a ( m ) J[ v e (z k )\ = 

n (* - n - %) l/ro n & - ^ m ■ ^ 

i,l k<l i<j 

In the next section, we present a comprehensive calcu- 
lation of the steps explained above for the non-Abelian 
FTSCs and identify the primary fields of the edge CFT, 



z = t + ix , z = t — ix , x ~ x + 2tt (6) 

where the periodicity of the x direction is due to the 
circular geometry of the edge (x is rescaled to satisfy 
the above relation and the Fermi velocity is chosen as 
vp = 1). By this transformation, the action of the edge 
theory becomes 

S = --^ J dzdz d z ip(z,z)d z tp(z 7 z) , (7) 

in which we have rescaled ip by a factor of \fm. We 
should be careful about this scaling as it has a physical 
effect. Before its action, the bosonic ip is assumed to be 
living on a circle with radius R = 1; i.e., <p is identified 
with ip + In. However, after it we should take R = y/m, 
i.e., ip ~ ip + Technically speaking the bosonic 

field is compactified on a circle with radius R = \fm. 

The equation of the motion for the above action is 
d z dz<p (z, z) — 0, so either d z ip (z, z) = or d z p (z, z) — 
0. We choose the latter choice, which means <p (z, z) — 
ip (z) since we need a chiral (holomorphic) edge state (the 
first choice given an antichiral (nonholomorphic) edge 
state. The edge state is enjoying a U(l) symmetry under 
which ip —></? + e with any arbitrary e. This symmetry 
leads to a conserved Noether current. It is given by the 
following formula: 

j (z) = id z p (z) , (8) 

which is manifestly conserved since d z j (z) = 
id z d z ip (z) = 0. The propagator of the ip field sat- 
isfies the following equation: 
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-^-d z d z K(z,z,w,w) = 6^(z-w). (9) 

Using the ^d z (z — w)^ 1 identity for the delta function 
in complex plane, we have: 

K (z, z, w, w) = (if (z, z) if (w, w)) = — In \z — w\ 

= -ln(z-w)-ln(z-w). (10) 

For holomorphic solutions of the edge theory we have 
(tp (z) <p (w)) = — In (z — w). Now let us find the correla- 
tion function of the current operators. 

(j (z) j H) - -d z d w (p (z) ^ H) - -^—5 (ii) 

(z — w) 

We can read the scaling (conformal) dimension of current 
operator off above correlation means leading to hj = 1. 
An interesting set of operators known as vertex operators 
can be obtained from the bosonic ip chiral fields. We have 

V a {z)=e ta ^ z \ (12) 
It can be easily shown that 

(V a {z)Vf } {w)) = (z-w) a r (13) 

An important fact about these operators is that they sat- 
isfy (V a (z) V H) = e™ af) (Vf) (z) V a (w)), which gives 
anyonic statistics for a(3 ^ Z. Taking a = (3 we can 
compute the scaling dimension of V a , which happens to 
be 

h a = a 2 II. (14) 

The fact that tp is compactified on a circle with R = ^/m 
allows only discrete values for a. We have 



current operator j we can measure its electric charge, 
which happens to be e. Therefore, V e is actually an elec- 
tron operator of U(l) CFT compactified on a circle with 
radius R = \fm. Additionally, by choosing a = ^= and 
computing the correlation function between two of them, 
exchanging them counterclockwise we obtain 9 = ^ for 
their fractional statistics. The charge of this operator can 
also be measured by looking at its OPE with the current 
operator and turns out to be e/m. Accordingly, e 1 ^ 
represents the Laughlin quasiparticle (anyon) operator. 



C. Chiral algebra of the Laughlin-type states 

In this section we wish to study the operator content 
of the Abelian FQH/FTSC states, through the so-called 
chiral algebra. The conformal algebra always contains 
the Virasoro algebra, which is given by a central charge 
c. The edge theory of the Abelian FQH/FTSC states at 
v = 1/m is described by a CFT with c = 1. The Viraso 
algebra at c = 1 allows for an infinite number of primary 
fields of the form V a j n with arbitrary real a and integer 
n numbers, where V a — e taip and j = id z p. In the phys- 
ical systems, however, we always have a finite number 
of primary fields. One way to obtain a finite number of 
primary fields is to extend the Virasoro algebra. A par- 
ticular example is through the chiral algebra where we 
choose one of the primary fields as well as the current 
operator j as the generators of the algebra and find all 
the primary fields that have nonsingular OPE with them. 
The resulting operators will form a representation of the 
chiral algebra. For example, in the c = 1 we can take the 
electron operator V e = exp (i^pmp (z)), , and j with 
scaling dimensions m/2, m/2 and, 1, respectively, as the 
generating operators. The OPE of the current operator 
with vertex operators V a is always nonsingular. Now let 
us consider the OPE of V e with vertex operators. We 
have 



V a (z) = e ia ^ z ' = e ia {<i> ( - z )+ 2lT V™) 
k 

a = —= , k e Z. 



(15) 



Note that there is no constraint on the value of the k 
and it can be any integer number. Therefore, we still 
have an infinite number of primary fields. In the next 
section, we discuss the chiral algebra that puts further 
constraints on the value of k. Now let us consider the 
following operator: 



Ve(z) 



_ e -iy/m<p(z) 



(16) 



The scaling dimension of this operator is m/2 and ex- 
changing two of them yields a minus sign. Accordingly, 
V e (z) represents a fcrmion operator. By looking at the 
operator product expansion (OPE) between V e and the 



V e (Z) V a (W) 



1 



(z — w) y 



(17) 



In order to get a nonsingular OPE, ^Jma should be an 
integer. Therefore, every operator where k e Z m is 

also allowed in the chiral algebra. These operators create 
excitations in the v = 1/m Laughlin states. The fusion 
algebra of this chiral algebra is pretty simple: V k /^i x 
Vij^ = Vk+i mod m , which represents a Z m algebra. 

The chiral algebra of the free boson theory compact- 
ified on a circle with R = ^fra radius can be further 
expanded by choosing the bosonic operators e ±4 ( 2 V™v) 
and j as the generators of the chiral algebra. This alter- 
native choice turns out to be very useful in studying the 
operator content of the U(l)/Z 2 orbifold theory descrip- 
tion of the superconductors. The intuition is that the 
charge of the Cooper pair is 2 and they are considered 



7 



as the boundstate of two electrons. On the other hand, 
e *(V^«) - lim^ V e (z) V e (to). Following the proce- 
dure we used above, we find the following operators as 
the allowed primary fields: 



V k/2Vm (z) = e l ^ v{z) ,h k = £,ke Z im . (18) 

Form the above relation, it turns out that the previous 
chiral algebra is a sub-algebra of the new one. The oper- 
ators that we encountered before and that represent the 
Laughlin state correspond to k = 21 and I < m. The 
new primary fields with k — 21, where m < I < 2m, were 
identified with k — 2m primary fields in the previous case 
(Z m algebra). The primary fields with k — 21 + 1 are 
all odd under the ip — > tp + 27Ty / m symmetry transfor- 
mation of the original free chiral boson theory, so they 
can be removed form the theory at the end. The re- 
sulting theory can be easily shown, which is modular in- 
variant, thus physical. The fusion rules for this chiral 
algebra reduces to a Z^ m one. In general, for a chiral 
boson theory compactified on a circle with a rational ra- 
dius of ~ = 3j form where p and p' are coprime (i.e., 
they their greatest common divisor is one), the maxi- 
mal chiral algebra is generated by j, and e ±Iv/ ^ v , where 
N = pp' , whose corresponding fusion algebra reduces to 
Z 2N group. [13 EI] This chiral algebra of the c = 1 
free boson theory is referred to as An ■ As a sanity check, 
let us consider R = y/m theory that represents the edge 
state of a v — 1/m Laughlin state, ^- = ^. Since m is 
odd (m,2) = 1, then N — 2m. Accordingly, the maxi- 
mal chiral algebra is given by A 2m a-ud the fusion algebra 
reduces to Z^ m . 



D. Vortices and the twist operators 

Imagine a superconducting state with a number of 
vortices. Because superconductor is a condensate of 
Cooper pairs with charge q = 2e , it allows for the vor- 
tices with a magnetic flux half of the quantum flux i.e., 
<t>o = ?p = n/e in the unit where h—c—1. On the other 
hand, the Aharanov-Bohm effect tells us that electrons 
moving around a vortex will pick up a minus sign. In 
other words, having a vortex at point R requires the in- 
troduction of a branch cut in the superconductor such 
that electrons change sign by passing through it. This 
means a vortex can be imagined as a topological defect 
in the superconductor. Now let us formulate this is in 
a more elegant way. We wish to consider an operator 

a (^Rj referred to as the twist operator which creates a 

vortex at point R. If we take an electron and turn it 
around this twist operator, the electron operator's sign 
changes. Within the bulk-boundary CFT theory corre- 
spondence conjecture this means we need an operator 



whose OPE with electrons is of the following form: 

V e (z)a(w)= aH (19) 
(z — w) 

Going around the twist field a (w), the electron operator 
picks up a minus sign. Therefore, the twist operator can 
change the boundary condition on the electrons from pe- 
riodic to antiperiodic and vice versa along space and/or 
time directions. Such a twist operator is quite well stud- 
ied in the context of the Ising CFT and physically corre- 
sponds to the spin order and disorder fields. [301 EH] 

E. Edge theory of the non-Abelian FTSCs: U(l)/Z 2 
orbifold CFT 

As we discussed before, it is known that the maximal 
chiral algebra of a free boson compactified on a circle 
with ^- = ^ with (p,p') — 1 is given by An, where 
N = pp'.jH |H E] This algebra is generated by the 
bosonic operators: j, and e ±l ^^ i P, Accordingly, the 

primary fields allowed in this chiral algebra are e % , 
where k € Z 2 n, iu addition to the current operator 
j = id z ip. Therefore, the maximal chiral algebra of a 
free chiral bosons compactified on a circle with R = ^/m, 
that describes an Abelian FQH (or an Abelian FTSC) 
state at v = 1/m filling, is given by Aim- The free boson 
edge theory enjoys a discrete Z 2 symmetry which acts on 
the U(l) boson field in the following way: 

K : cp ->■ -cp, (20) 

or, equivalently, 1Z : j — > —j. The action of the 
Z 2 group 1Z divides the primary fields of the A 2m chi- 
ral algebra into the invariant and noninvariant groups: 
4>k = cos (ktp/2y^rn), and 4>k — sin (k(p/2y / m). On the 
other hand, as we discussed in the previous section, for 
a range of parameters adding superconductivity to the 
system can gap out half of the degrees of the freedom 
corresponding to primary fields of the edge theory and 
keep the rest gapless, meaning either the invariant or the 
noninvariant part of the primary fields should be taken. 
As a result, the electron operator can be taken as the 
cos (y/mip), which is even under the Z 2 orbifold theory. 
The above argument suggests that we should look 

for a chiral algebra generated by j, and cos (^/2N(j)^ 

(N = 2m). The resulting CFT is known as the U(l)/Z 2 
orbifold theory that is achieved by moding out the dis- 
crete symmetry 1Z from the free boson edge CFT. More- 
over, the chiral algebra of this orbifold theory is given by 
A N /Z 2 . 

The orbifold theory is known to be relevant in other 
topological states of matter as well. In Rcfs. {31 12"51 IM] . 
Brakeshli and Wen studied a bilayer FQH state with a 
strong interlayer repulsion. In the absence of strong re- 
pulsion, the edge theory of the bilayer FQH state is given 
through two free chiral boson fields ipi = Lp+J ^ p ~ and 
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<~pi = ip+ t, ip ~ and thus enjoys a £7(1) x £7(1) symmetry 
corresponding to j + = id z ip+ and ?_ = id z (p~ conserva- 
tion. A discrete symmetry of the edge theory corresponds 
to the exchange of the two layers: 1Z : if- — > —<p~- In the 
(to, to, 0) bilayer FQH, the electron operator in each layer 
can be represented as ipi = e l V™v>i ^ where i = 1,2. By 
increasing the interlayer repulsion, it is energetically fa- 
vorable for the quasiparticles to reside in the same layer. 
This causes only the even combination of electron opera- 
tors under the 1Z symmetry to remain gapless and the odd 
combination becomes gapped. Consequently, the edge 
CFT is generated by = e^™^/ 2 cos (y/rrvp-/2). 

Therefore, the edge theory is given by the free £7(1) the- 
ory, corresponding to ip+, and a U(l)/Z 2 for the the im- 
part, thus a £7(1) x U{l)/Z 2 orbifold CFT. As in our ex- 
ample, the orbifold theory allows for the twist fields to 
show up, resulting in a non-Abelian state. In a different 
study, Barkeshli and Qi have shown that Chern insu- 
lators with higher Chern numbers can provide a physi- 
cal playground to realize a similar physics. [38] More re- 
cently, Barkeshli et al. have studied the topological order 
and projective non-Abelian properties of these types of 
models. 1 32 



F. Operator content of the U(l)/Z 2 orbifold CFT 

To obtain the chiral algebra representing the non- 
Abelian FTSCs, we should mod out the Z 2 subgroup of 
the CFT that acts on the current operator as 1Z : j —> 
—j. Moding the Z 2 symmetry has two effects. First, 
the resulting primary fields should be either symmetric 
or antisymmetric under the orbifold action. Second, it 
allows new primary fields to show up in the chiral alge- 
bra. These new primary operators are twist fields. The 
OPE of any primary field in chiral algebra with the gen- 
erators of the chiral algebra, i.e., j and cos (^/2Nip S j , 
should be nonsingular. The OPE of the twist fields with 
the cos (y / 2Ntpj vertex operator is nonsingular. How- 
ever, its OPE with the current operator j = id z ip has 
a branch cut, which means taking the current operator 
around the twist field changes the sign of the current op- 
erator leading to a singular behavior of the OPE. For the 
unorbifold theory with the full £7(1) symmetry, this is 
not acceptable and as a result twist fields are absent in 
the chiral algebra. However, by moding the Z 2 symme- 
try generated by 1Z, the twist fields are allowed because 
—j (z) = TZ j (z) TZ ~ j (z). In other words, the current 
operator is invariant after going around the twist fields 
up to the action of the Z 2 group. Therefore, the chiral 
algebra of the An/Z 2 theory has to include twist fields 
as well. It is these twist fields that make the orbifold 
theories significantly richer and more exciting than the 
original An chiral algebra associated with the Abelian 
FTSC. 

There are two kinds of twist fields in the £7(1) jZ 2 orb- 
ifold theories: o\ and a 2 fields with conformal dimension 



h a = 1/16 and n and t 2 fields with conformal dimen- 
sion h T = 1/16 + 1/2 = 9/16. The latter twist fields are 
related to the first one through the following relation: 



j (z) <7i (w) 



(z — w) 



(21) 



which is consistent with the fact that cr^ (w) twist opera- 
tors change the sign of the j (z) — id z (p (z) field through 
27r rotation around w. As was stated before, twist opera- 
tors are allowed by the chiral algebra An /Z 2 since under 
the action of the twist field on the current operator, it is 
invariant up to the action of the Z 2 group (more precisely 
j = — j in the Z 2 orbifold theory). 

The generators of the chiral algebra are the current (j) 

and cos ^y / 2Nip S j operators (N = 2m). We now consider 

the rest of primary fields that have nonsingular OPE with 
these generators other than twist fields. It can be shown 
that the following set of operators has nonsingular OPE 
with these operators. 



(22) 



which are simply the ^-invariant part of the primary 
fields of the A 2m chiral algebra that we encountered be- 
fore. There are two fermion operators in the A 2m /Z 2 
chiral algebra among which we take the ^-invariant one 
as the electron operator. They are 



JV 




sin [y/m(p 



(23) 



Several properties of these primary fields are summa- 
rized in table 



III F| . The most important information 



one needs in any CFT is the fusion algebra of the primary 
fields. For the Z 2 orbifold theories we are concerned with 
here, Dijgkgraaf et al. [22) have derived their fusion rules. 
According to their study, for even (in our problem it 
is the case) the fusion algebra is as follows: 



4>k x 4> k ' = <pk+k' + <Pk-k' (k' 7^ k,N - k) 
4>k x 4>k = i + j + (t>2k 
4>N-k x 4>k : 



<hk + 4>N + <fir- ( 24 ) 



The above fusion rules can be verified by interpreting 
4>k fields as cos (k(p/ {^/2N^j. Also, the fusion rules 



involving j, 



b N and 



operators are: 



j x 4> k = <fik, 
jxj = 1, 



4> N x ft = 1 



!J N 



J N 



3, 



(25) 
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TABLE I: Primary fields of the orbifold theory given by an 
An/Z2 chiral algebra with N = 2m, and the corresponding 
conformal (scaling) and quantum dimensions. 





CFT primary 


Conformal 


Quamtum 




field 


dimension 


dimension 





1 





1 


1 


3 


1 


1 


2 


4>2m 


to/2 


1 


3 


<fik 


k 2 /(8m) 


2 


4 


(Ti 


1/16 


a/2to 


5 


Ti 


9/16 


V2m 



important observation in studying the zero modes of the 
anyon operators. 

IV. TOPOLOGICAL ORDER OF THE 
NON-ABELIAN FTSCS 

In this section, we study different properties of the 
non-Abelian excitations. We prove the existence of the 
zero modes bound to vortices, compute their quantum 
dimensions and derive the non-Abelian part of the braid 
matrix that acts on the degenerate ground state by braid- 
ing fractionalized Majorana zero modes. 



and, most importantly, the fusion rules involving the 
twist operators is given by the following formulas: 

(Ti X (Ti = 1 + (f) l N + ^ <fik, 
k even 

0-1X0-2= ^ <f>k, 

k odd 
j X (Ti = Ti, 

(fik X (Ti = <7j, 

4> l N x a i = (26) 

The fusion rules for the twist fields t$ can also be read 
off the Ti = j x ai formula. The above relations suggest 
that the AzmjZi chiral algebra supports 2m + 7 primary 
fields, among which (jy^, (fi 2 N , and j are Abelian fields 
and (frk's are non-Abelian primary fields with a quantum 
dimensions equal to 2 and the non-Abelian twist fields 
with a quantum dimensions equal to %/iV. Now let us 
make comments on the above fusion algebra. First of all, 
as we stated before, (fi\ m field with conformal dimension 
h = to/2 is a fermion field and invariant under under the 
1Z symmetry. We choose it to represent the electron oper- 
ator. The above fusion rule between an electron operator 
and a twist field a± along with their conformal dimensions 
suggest the following OPE (and a similar relation for the 
ri): 

V & (z)a(w) = ° {W) (27) 
[z — W) 

This relation is consistent with the fact that taking an 
electron around the twist field flips the sign of the elec- 
tron operator and changes the boundary conditions. Sec- 
ond, (fin — cos (^J^j operators represent the 1Z invariant 

non-Abelian anyon operators with conformal dimensions 
h 2 i = I 2 /2m. The OPE of these anyon fields with the 
twist fields o\ is (and a similar relation for the t\) 

(fi 2l {z)a{ W )= a( g (28) 
(z — w) 1 

which means taking an anyon operator around the twist 
field will change its phase by 621 — Pn/m. We use this 



A. Mode expansion and zero modes 

In this section, we study the physical effects of the twist 
fields and show that they allow for the existence of the 
fractionalized Majorana zero modes. First of all, let us 
mention the following mode expansion for the primary 
fields. For an arbitrary chiral primary field $ (z) with 
conformal dimension /i$ we can write the following mode 
expansion: 

*(*) = E ^ n',p,qGZ. (29) 

n—n'— — 

Q 

In order to determine the allowed values for p and q 
integers, we need to know the behavior of the $ (z) field 
under the twist, i.e., the value of the 9 in $ (e 27Tl z) — 
e $ (z). On the other hand, the above mode expansion 
suggests that $ (e 27ri z) = e -«(^-p/?) 2 ^$ ( z ). Therefore, 
I = /i$ — ^ mod 1. If p/q = 0, then it means cto 
is allowed in the system and we say that the & field 
supports a zero mode. For physical systems we usually 
have periodic boundary conditions on fields, i.e., 8 = 0, 
and therefore for the primary fields with rational con- 
formal dimension zero mode is absent. For example, 
in the Laughlin FQH state at v = 1/m filling fraction, 

q = f - h mod 1 = V 2 - This means n = m + 1/2, 
and accordingly, there is no electron zero mode. An in- 
teresting question that arises is what can cause a twisted 
boundary condition (i.e., 8 7^ 0) on the $ (z)l The an- 
swer is simply the twist fields. A physical realization of 
the twist field is the flux insertion (vortex) that through 
the Aharanov-Bohm effect can change the boundary con- 
dition. For example, a single half vortex with (fio/2 mag- 
netic flux can change the boundary condition on electrons 
from a periodic to an antiperiodic one. For such a bound- 
ary condition, | = ~ — ^ mod 1=0. Consequently, 
the zero mode is now allowed. This is, for example, the 
reason why zero modes are bound to the vortices in p+ip 
superconductors. 

As the fusion rules of the U(X)jZi orbifold theory sug- 
gest (see Eq. [26]), fusing two twist fields (vortices) of 
the same kind can generate 02fe primary fields. Now let 
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us take the anyon primary field with the least conformal 

dimension i.e. 02 = cos (^y=fj among them. The mode 

expansion of this primary field whose conformal dimen- 
1 

2m 



sion is 7^- reads as follows: 



n—n' — 



^n+l/2m 



n',p,qe z. 



(30) 



In the usual periodic boundary condition, p/q = 1/m, 
and therefore n = n'—l/(2m), so the zero-mode operator, 
ao, is absent. Now imagine we have put a twist field 
corresponding to the half vortices at w = 0. Taking the 
4>2 (z) field around the twist field will result in 6 = ^. As 
a result, in this case | = ^ — ^ m °d 0. Accordingly, 
n G Z and the zero-mode operator ao exists in the mode 
expansion of the 02 (z) field in the presence of a twist 
field at w — (or by flux insertion at w — 0). 

Let us consider the fusion algebra between <f>k and 
<j>2N-k fields, which is 0^ x (f>2N-k = 1- This means 
4>2N-k is the conjugate of the <pk] i.e., their fusion gives 
the vacuum. Thus, we represent the <j)2N-k as <jy k . In 
particular, take the anyon operator, <p 2 , whose conju- 
gate field is 02(m-i) w ith a conformal dimension equal 

to ^ 2? ^~ 1 ' > ■ According to their fusion, the OPE between 
these fields is given by the following formula (up to a 
constant number): 



4>2 i Z ) ^2 M 



1 



(z — w) 



2(m-l)+l/m : 



02 H 4 (z) ~ 2(m _ 1) + 1/m , 

(w — z) y ' ' 



(31) 



from which we obtain the following commutation rela- 
tion: 

4 (z) 2 (w) = e-™' m 4>2 M 4 {*) • (32) 

Using the mode expansion, the commutation relation 
between the modes can be achieved, which is 



(33) 



Another important property of the 02 field is that fusing 
it 2m times with itself gives the identity, i.e., 0| m ~ 1. 
This property has an important consequence for the zero- 
modes. To see that we first mention that modes can be 
written in terms of the following contour integrals: 



n^f^M*)^ 1 '™- 1 , (34) 

where c encloses the origin. Using this relation, it can be 
shown that 



The above relation suggests that ao is the generalization 
of the Majorana zero ipo mode for which ipfi = 1. In this 
problem a™ = aj m acts like a Majorana zero mode and 
for mo > 1 we refer to the zero mode operator «o as the 
fractionalized Majorana zero mode for obvious reasons. 



Quantum dimension of the non-Abelian 
excitations 



One of the features of the topological order is the GSD. 
This degeneracy is related to the so-called quantum di- 
mension of the zero modes. A zero mode, as its name 
implies, excites quasiparticles with zero energy. In other 
words, by applying a zero-mode operator on the ground 
state we yield a new state with the same energy. The de- 
generacy of the ground state is a topological index of the 
system which is robust against perturbations. Clearly, 
the ground state(s) of a Hamiltonian depends on the 
boundary conditions. The existence of some fields can 
change the boundary condition and as a result the GSD. 
For example, as we mentioned in detail before, a twist op- 
erator at the R\ point exchanges the boundary condition 
on electrons moving around it and gives the possibility of 
the zero modes. One way to define the quantum dimen- 
sion of an operator, in this case twist fields, is to consider 

n of such operators at isolated points ^R\, Rr^j ■ Then 

compute the GSD with respect to this constraint. The 
GSD usually grows with the number of such fields. Ac- 
cordingly, the quantum dimension of these operators is 
defined as follows: 



lnd 



\nGSD 



n 



CO. 



(36) 



a 



2m 



J 277 



1. 



(35) 



There are several other ways to define the quantum 
dimension of an operator, but we stick to the above defi- 
nition for our purposes in this paper. Now let us consider 
an FTSC on a torus (or an infinite plane) with two vor- 
tices at R\ and i? 2 in it. The vortices can be modeled 
by putting a\ twist fields at these points. As we argued 
before, twist fields allow the existence of zero mode oper- 
ators bound to (localized at) twist operators. This bound 
state can be justified in an intuitive way. Naively, a vor- 
tex means a circular puncture with radius L (of the order 
of the correlation length of Cooper pairs) centered at Ri . 
We expect an edge state propagating along the boundary 
of the puncture and localized in the perpendicular direc- 
tion. The reason for this expectation is that inside the 
puncture is a trivial state while outside it is an FTSC. 
The chiral anyons (corresponding to the 02 field ) that 
propagate along the boundary have energies of the form 
E = v F k, where k is their momentum and w F represents 
the Fermi velocity. However, since L is finite, the mo- 
mentum of the quasiparticles is quantized as k n = ^j^, 
where n is an integer due to the fact that the twist field 
(vortex) allows for the zero mode. Therefore, level spac- 
ing is AE oc l/L — > oo, meaning that only the zero mode 
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(n = 0) can be observed at that point. Therefore, twist 
fields at Ri allows for the fractionalized Majorana zero 
mode, ato, at that point and similarly does the twist field 
at i?2- Let 71 and 72 operators denote zero modes at R\ 
and R.2- They follow the following algebra: 

7172= exp H7r/m) 72 7i , 7?™ = 7 2 2m = 1(37) 

Using the above algebra we can show that 7^72 operator 
satisfies the relation 

/ , \ 2m 

(7I72) =exp(-*7r(2m+l)) = -l, (38) 
which defines a 2m state with the following eigenvalues: 
e -/(2m) 7 t 72 \ q) = exp faq/m) \ q ) . (39) 

Therefore, we can consider e ln '/( 2m ) 7^72 as the general- 
ization of the familiar fermion parity operator for Majo- 
rana fcrmions. Hence, we call it the i?2m character of the 
ground state. To understand the nature of the degen- 
erate ground states, the following relations that can be 
easily proven will be found useful: 

7} \q) = exp (-«r£i (q - 1/2) /m) \q + 1 mod 2m) , 

71 \q) = exp (in^i (q — 1/2) /m) \q — 1 mod 2m) , 

72 \l) ~ ex P (~* 7r C2 (q — 1/2) /m) |<7 + 1 mod 2m) , 
72 \q) = exp (i< 2 (9 - 1/2) /m) |g - 1 mod 2m) (40) 

where £2 — £1 = 1, and £2 can be any real number. The 
ambiguity in the value of the free parameter £ 2 indicates 
that only quadratic combinations of the form 7^72 are un- 
ambiguous, and thus physical and measurable. However, 
the value of the £2 does not affect the quantum dimension 
of the fractionalized Majorana zero modes (7, operators) 
as well as their projective braid statistics. Therefore, we 
fix the gauge and assume £2 = 1> and £1 = as of now. 
With this gauge convention, the ground state \q) can be 
constructed in either of the following ways: 

\q) = 7l 9 |0) = exp (m/2m (q 2 - 2q)) |0) . (41) 

The above relations play a key role in our paper. To see 
how powerful the CFT methods is, note the fusion rule 
given for two twist fields in Eq. [26] - It simply means 
that taking every two twist fields (vortices), their com- 
posite state results in different fusion channels which is 
a signature of their non-Abeilan nature. The final state 
can be either vacuum or a state with one of the anyon 
operators 02, <f>±, and <f>i m -2- Note that the zero 
mode of the 4>2k operator can be created by taking the 
&-th power of the fractionalized Majorana zero mode of 
the c/>2 i-e. 7 2 fc (i can be cither 1 or 2 since we have fused 
twist fields at R\ and R2 ) , which can be easily proven us- 
ing the fusion rules for 2 z operators. Therefore, fusing 
every two twist operators results in one of the following 



states j|0) ,tJ |0) , 7 | 2 |0) , ...^t 2 ™" 1 |0)|. As a conclud- 
ing remark, remember that fractionalized Majorana zero- 
mode operator 7$ is bound to o\ (Ri). Therefore, fusing 
two twist fields can be viewed as considering a 7^72 op- 
erator, which happens to define a 2m degenerate states. 
So the results of the CFT are consistent with the direct 
calculation. 

The main result of our discussion so far is that by tak- 
ing any two vortices, we can define the Z2 m character 
e nr/(2m)^\^, ^Jiat takes values in Z2 m and as a result 
defines 2m degenerate states. Accordingly, taking n vor- 
tices, we can define n/2 of such ^ m characters (remem- 
ber that in a compact space there are always an even 
number of vortices) and results in the (2m) n ^ 2 degen- 
erate states. However, we should note that the num- 
ber of electrons in a superconductor is conserved mod 2. 
Each electron operator is formed of m anyon operators 
(4>2 fields) in A2m/Z 2 chiral algebra. So the number of 
anyons is conserved mod 2m. This generalized parity of 
anyons can be computed through the characters. It 

is indeed x — Y\i=i el7r ^ 2 "^7;L-i72i- The conservation of 
the x operator reduces the number of degenerate states 
by a factor of l/(2m), since one of the Z 2m characters is 
fixed by the value of the x and the rest of i?2m characters. 
Altogether, we reach the following conclusion: 

GSD (n) = (2m) n/2_1 => d a = V2m. (42) 
which agrees with the result of the CFT. 

C. Projective braiding statistics of the 
fractionalized Majorana zero modes 

Braiding means an adiabatic evolution of the ground 
state in space-time, after which two quasiparticles ex- 
change their position. For example, exchanging the po- 
sition of two electrons causes the change of the ground 
state wave function by tt. In the FQH state at v = 1/m, 
exchanging the position of two anyons counterclockwise 
changes the phase of the ground state wave function by 
7r/m phase factor. For non-Abelian quasiparticles, how- 
ever, the ground state is degenerate and exchanging two 
quasiparticles may cause a unitary transformation of the 
ground state; i.e., it may mix different degenerate ground 
states after exchanging two non-Abelian quasiparticles 
(for example two Majorana zero modes) adiabatically. In 
the CFT, braid matrix (that acts on the ground states) 
is defined through the four-point correlation functions. 
Take four primary fields $ (21), $ (z 2 ), $ (z 3 ), and $ (24). 
Then compute the following correlation function: 

($ (zi) $ (z 2 ) $ (z 3 ) $ (z 4 )) • (43) 

Suppose the fusion algebra of the $ primary fields is given 
by $ x $ = 53i =1 ST;; for some know fi; fields. Using 
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this fusion algebra, there are several ways to compute 
the four-point correlation function. One way is to fuse 
$1 with $2 and $3 with $4 in the I channel and sum 
over different channels. For each particular channel let 
us call the result Ffn (l,z), where z — Zl2Z34 and za — 
Zi~Zj. The second way is to fuse $1 with $3 and $2 with 
$4 in the V channel and sum over different V channels. 
Again for each particular channel let us call the result 
JF\ 3 (l 1 , 1/z). The consistency check requires T\ 2 z ) to 
be a linear combination of T\ 3 (I 1 , 1/z) functions say: 



a, = 



23j; 



(44) 



-B23 is called the braid matrix associated with braiding 
$2 and $3, which acts like a unitary transformation on 
the degenerate ground states formed of four $ field in- 
sertions. Let us go back to our problem. Upon vortex 
(the dual of the edge CFT's twist fields in the bulk of the 
FTSC) insertion, a fractionalized Majorana zero mode is 
created at that point. Instead of braiding zero modes, 
we can instead braid vortices that obviously braid the 
zero modes bound to them as well. It should be em- 
phasized that the braid matrix can be decomposed in 
the Abelian (phase factor) and non-Abelian parts. The 
Abelian part depends on the microscopic details of the 
system and the interaction between superconducting vor- 
tices. Since it does not affect the non-Abelian part of 
the braid matrix, pQ we neglect the Abelian part of the 
braid statistics in this paper and refer to the remaining 
non-Abelian part as the projective braid matrix. [32] The 
calculation of the projective braid matrix from the con- 
formal blocks using the above recipe is straightforward 
but tedious. In this section, we employ a somewhat eas- 
ier but related method that uses the S matrix to compute 
the non-Abelian part of the -B23 matrix. To this end, we 
first compute -B12, which acts on the ground state upon 
exchanging twist fields <j\ and <r 2 (as well as the zero 
modes bound to them) . To compute -B12 we can actually 
forget about other twist fields and consider vortices 1 and 
2 only. From the fusion rule given for two twist fields in 
Eq. [26] , we can write the OPE between a\ and a 2 fields 
(up to constant coefficients). We have: 



(71 (z) (7l (w) 
1 

(z — w) 



2m- 1 

(z-w) 

L 9=0 



q 2 /(2r, 



W + Z 



.(45) 



exchanging z and w [accordingly a\ (z) and a\ (w)] in 
the qth channel, one can verify from the above OPE that 
it results in an e l7r / s e ~ l7: <i / 2m phase factor. Forgetting 
the Abelian common factor in every fusion channel, i.e., 
e"/ 8 , and keeping the g-dependent part only, we con- 
clude that in the q-th channel, i.e., \q) — j\ q |0) exchang- 
ing 71 and 72 as a result of exchanging twist fields, results 



m e 



-iirq /2m 



\q). Therefore, 



Bia |g) - [B 12l \B\ 2 ) B 12 |0> = 

exp (-mq 2 /2m) \q) = exp (-iwq/m) |0) , (46) 

stating that the -Bi 2 operator is diagonal in this basis. 
From Eq. 
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, we can rewrite the above formula as: 
B 12 \q) = exp {-mq/m) 7 2 t<? |0) . (47) 



Since B 12 |0) = |0), by comparing Eqs. [IT] and [47] 
we find that B 12 ^ x b\ 2 = e"/ m ^ 2 - To find the way the 
braid operator acts on 72 we can use a similar procedure. 
We can also use the fact that the Z 2m character, i.e., 
the e I7r ^ 2m ^~f\~f 2 operator, is a physical and measurable 
quantity and cannot change after braiding 71 and 72 and 
keeping the rest of the system untouched. To guarantee 
this physical expectation, the braid operator has to act 
on 71 and 72 as follows: 



^1271^12 



— p i7r / 7 



B 12 j 2 B 



t _ p i7z/m 



72, 
t. 



727i72- 



(48) 



Now we wish to compute the -B23 braid matrix from the 
diagonal -B12 braid matrix. For this purpose, we have to 
consider four vortices at least. Imagine we are given four 
vortices at Ri points, where i = {1, ..,4}. Taking the fist 
two vortices, we can define a |g) 1 2 state that diagonalizes 

the e 17 * '/( 2m ' 7^72 operator as well as B12. We can do a 
similar job for the third and fourth vortices and define 
a \q') M state on which the e m ^ 2rn ^\^^ and S12 opera- 
tors act diagonally. Therefore, the ground state of the 



whole system is \q, </) 



(1,2)(3,4)' 



However, we should em- 



phasize that q and q' are not independent as the total Z2 m 
character of the system xt = e l7r ^ 2m ^7|72e i7r ^ 2m ^7374 = 

giTr 2 ^- is fixed. For simplicity, let us assume xt = 1, so 
q + q' = mod 2m. So we have: 



Remember that 4> 2m = (f> 2m = V e and <fio = 04 m = 1. 
The above formula means that taking two twist fields 
and fusing them (and as a result fractionalized Majorana 
zero modes bound to them) is ambiguous and depends 
on the fusion channel. The qth fusion channel contains 
k fractionalized Majorana zero modes. Therefore, ex- 
changing two twist fields does not simply give us a phase 
change like Abelian quasiparticles. Instead, the phase 
factor that we obtain depends on the fusion channel. By 



B i2 W, 9)(i, 2 )(3,4) = ex P (-inq 2 /2rn) \q, g) (li2 )( 3 ,4) . 

B34 \q, 9)(i, 2 )(3,4) = ex P (-W7 2 ™) k, <7>(1,2)(3,4)( 49 ) 



To find B23 w e need a basis in which e J7r ^ 2m ^7372 and 
as result -B23 is diagonal. Since xt = 1 is conserved, 
e i7r/(2m)^,t^, 4 w jjj k e diagonal on that basis as well. In 
other words, we need states of the form \p,p), 2 ^ in 
which B23 acts as e~ l7Tp l 2m . So the procedure is simple. 
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Through basis rotation, we can diagonalize B23 and then 
we rotate back to the former bases. Let us assume there 
is an S matrix that relates the two bases. This is actually 
the same S matrix that we encounter in the CFT on a 
2D torus that exchanges the two nontrivial cycles of the 
torus. Let us find the S matrix by simple calculations. 
We have 



|P'P)(3,2)(1,4) =J2 S p\ P >P\h 



2)(3,4) 



(50) 



The above relation suggests that e'^^ 27 ™^^ acts diag- 
onally on the \p,p)/ 3 2 )(i 4) state with an eigenvalue equal 
to e " p/m . On the other hand, 

in /(2m) t \ n tA _ 
e 7372 I"7j<Z/(3,2)(1,4) — 

y^^p \P~ 1;P+ 1)(1,2)(3,4) = X] ^P+l l^'^(l,2)(3,4) ' 



(51) 



To have a consistent result, we should have: S^ +1 
exp {iqir/m) which can be solved as follows: 



exp (ipqir/m) 



(52) 



To go from the |p,f>) (3;2)(M) to the k, 9>(i, 2 )(3,4) basis, 
we can use the >S _1 matrix. It can be easily verified that 
S 1-1 = 5*. Now let use start from the \q, q)^ 2 ^ 3 ^ state 

and braid vortices 2 and 3. To this end, we fist use to 
go to the |p,p)( 3 2 )(i 4) basis. Then we act on this state 
by B23, which is diagonal now. We finally rotate back 
to the 2 )(3 4) basis. These steps can be seen more 

explicitly in the following relations: 



•823 k,<?)(i i2 



)(3,4) 



VS^B 23 \P,P) 



(3,2)(1,4) 



Sl p exp(-ip 2 ir/2m) |p,p) (3i2)(M) 



q.k 

^(S 2 3^|fc,fc) (12)( 3 i4) , 



which can be summarized as: 
B23 = S^B 12 S , (Bi 2 )l 



-iq 2 Tr/ (2m) 



(53) 



Since B23 is not diagonal in the \q,q)n 2 )(3 4) basis, 
it means that exchanging vortices 2 and 3 (or 1 and 
4) causes a unitary non-Abelian transformation on the 
ground state. In other words, it will creat other degen- 
erate ground states upon braiding vortices 2 and 3. In 
contrast, after braiding vortices 1 and 2 (or 3 and 4), the 
ground state just acquires an Abelian phase which de- 
pends on the ground state though. An interesting ques- 
tion that arises concerns how braiding can change the 



ground state despite its adiabatic nature. The answer 
is that all the degenerate ground states share exactly 
the same physically measurable quantities. For example, 
they all have the same Z 2m characters and the same en- 
ergies. Additionally, having vortices in the system causes 
branch cuts. Upon braiding, some of the zero modes that 
bound to vortices have to cross these branch cuts. As a 
result, the zero modes at two vortices involved in this 
crossing will rearrange. This rearrangement modifies the 
wave function of the many-body system and the result- 
ing wave function should be still be in the ground space, 
meaning that it is a linear combination of different de- 
generate sectors of the ground space. 



V. WAVE FUNCTION OF FTSCS 

As we discussed in Sec. Ill, according to the bulk- 
boundary correspondence conjecture for the chiral states, 
e.g., FQH states, the wave function of the bulk Hamil- 
tonian can be written in terms of the conformal blocks 
of the edge-state CFT in an appropriate way. The first 
step is to identify the CFT that describes the edge of 
a chiral state. For an FTSC whose parent (unpaired) 
state is a Laughlin state at v = 1/m, as argued in Sec. 
Ill, the edge CFT is given by the A2 m /Z 2 chiral alge- 
bra of a U(l)/Z2 orbifold theory. The primary fields of 
this CFT contain two fermion fields <j)\ m = cos (y/rruf) (z)) 
and <j)\ m — sin (y/m<f) (z)). These are the counterpart of 
the Majroana fermion fields for the m — 1 case. In the 
FTSC, one of these two operators become gapped and 
the other remains gapless. Let us assume that cf)\ m is 
gapless. The wave function in the absence of vortices for 
a fixed number of electrons with positions at {zi}, where 
i = {1, .., N}, can be written in the following way: 



N 



* ({*}) - (ti m (*l) -41m = Pf (tTT^ 



Zj) 

(54) 

The BCS-like wave function can be easily obtained by 
summing over a different number of electrons (with the 
same parity as that of N). The result would be: 



\g) = exp 



(Zi - Zj) 



1 ,v^~ 



|0>. 



(55) 



Now let us find the wave function in the presence of 
2n vortices. Because the twist operator a\ {r\) inserts a 
vortex at rj, the wave function for n vortices at rjj where 
i = {1, .., 2n} and for a fixed number of electrons at {z{\ 
points where i = {1, ..,N}, can be similarly written as: 



* {{vj} ; {*}) = (<n (Vi) -<?i (Vn) <f>l m (zi) -4\ m (zn)) ■ 

(56) 
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The above correlation function can be obtained by us- 
ing the OPE between twist fields. As Eq. 
two twist fields can be fused in 2m different channels 



states, 
In 



other words, although the positions of the vortices are 
known, the above correlation function can be evaluated 
in different ways. We can divide the n vortices in n/2 
pairs, fuse every two twist fields in a pair, and then find 
the correlation function. Again, we should be careful 
about the last pair, because after all fusions, we should 
end up with identity operator for a non-vanishing corre- 
lation function. Therefore, the last channel is dictated by 
the rest of the system. Therefore, there are (2m)"^ 2_1 
different conformal blocks, each corresponding to a wave 
function corresponding to the {<?i, q n /2-i } choice 

for fusion channels. Therefore, the number of allowed 



degenerate ground states, i.e., the GSD, is (2m 
n vortices in the FTSC. 



.n/2-1 



for 



VI. Z 2m GAUGE THEORY OF THE 
FRACTIONAL TOPOLOGICAL 
SUPERCONDUCTORS AT 1/m FILLING 
FRACTION 

Here, we analyze the gauge symmetry of the v = 1 /m 
FTSCs and make a connection between non-Abelian 
FTSCs and the Zm rotor model discussed in Ref. [35]. In 
the FQH state, because of the number conservation, the 
gauge theory is described by a U(l) gauge field. However, 
in the presence of pairing, this symmetry breaks down, 
and the number of electrons is conserved only modulus 
2. On the other hand, the electron operator can be imag- 
ined as the bound-state of m anyon operators: V e = 4>™. 
Therefore, adding or removing two electrons is equivalent 
to adding or removing 2m anyons. In other words, the 
total number of anyons is conserved only modulus 2m. 
Therefore, if we imagine anyons as the fundamental exci- 
tations of the ground state, the low-energy theory of the 
FTSC at v = 1/m filling fraction is invariant under the 
Zi m gauge transformation. On the other hand, as we dis- 
cussed before, the quantum dimension of the topological 
defects (vortices) is \/2m, which agrees with the predic- 
tion of Z 2m rotor model. Using the latter formalism, You 
and Wen have shown the quantum dimension of the dis- 
locations (vortices) in the Zjq model to be y/N. [35] From 
this consistency check, we speculate that Zi m symmet- 
ric FTSCs at 1/m filling fraction can be thought of as a 
realization of Zi m rotor model. 



VII. SUMMARY, CONCLUSION, AND 
OUTLOOK 

In this paper, we introduced the notion of an FTSC. 
There are two classes of FTSCs, Abelian and non- 
Abelian. The Abelian FTSC is topologically the same 
as the Abelian FQH state, except that the electron num- 
ber conservation is violated by injecting/inducing Cooper 



pairs into the FQH state. The superconducting order 
parameter is small in the Abelian FTSC so that we can 
smoothly evolve the Hamiltonian of an Abelian FQH into 
that of an FTSC without closing the energy gap. How- 
ever, the non-Abelian FTSC is achieved upon closing 
and reopening the many-body energy gap of an Abelian 
FTSC by increasing the pairing amplitude. The result- 
ing state allows non-Abelian FMFs with quantum di- 
mensions equal to \j2m, where 1/m refers to the fill- 
ing fraction of the parent FQH state. We studied the 
braid statistics of these non-Abelions and derived the as- 
sociated projective braid matrix. The wave function for 
the non-Abelian FTSCs was suggested through comput- 
ing the conformal blocks of the edge CFT. We finally 
discussed the gauge theory of the effective low-energy 
Hamiltonian of the non-Abelian FTSC which happens to 
be a Zi m one. This observation makes a clear connec- 
tion between our model of the non-Abelian FTSC with 
the Zim rotor model with similar topological properties. 

Due to the robustness of the zero modes in the topo- 
logical phases against local perturbations of the Hamil- 
tonian and disorder, the non-Abelian anyons are promis- 
ing candidates for the topological quantum computation. 
Information can be stored nonlocally in the degenerate 
ground states, and braiding of quasiparticles causes uni- 
tary transformation on the ground state so can be used 
as the quantum gates. However, it has been shown that 
when the square of the total quantum dimension of a non- 
Abelian system is integer, it cannot perform universal 
quantum computation by braiding only, [31)] unless sup- 
plemented by measurement at the intermediate stages 
of the phase gate. [31 [37J For non-Abelian FTSCs orig- 
inated form Laughlin state at v = 1/m, the square of 
the total quantum dimension, D 2 = J2i=i d 2 , where N qp 
is the number of quasiparticle, and di their quantum di- 
mension, is integer. Hence, they are capable of universal 
quantum computation with braiding and measurement. 

It is worth mentioning that the heat conductance in 
the FTSCs depends only on the difference between the 
central charge of the right- and left-mover boundary the- 



ories. Accordingly, we have k x 



T . This equa- 



xy ~~ 67rfi _ 

tion can be proven by deriving the gravitational Chern 
Simons action of the FTSCs. 

Finally, we would like to comment on the topological 
field theory (TFT) of the bulk of a non-Abelian FTSC. 
Moore and Seiberg [32] have shown that for a system 
whose edge theory is given by G/H orbifold CFT, where 
H is a discrete subgroup of G, its bulk TFT is a Chern- 
Simons (CS) theory with G x H gauge group. In our 
case, G = U(l) and H = Z 2 , so the bulk TFT of a non- 
Abelian FTSC at v = 1/m is described by an 0(2) CS 
action at level m. Furthermore, the gauge symmetry of 
the v = 1/m FTSCs was discussed in this paper and 
shown to be a Zi m gauge symmetry. This observation 
makes our prediction about the quantum dimension of 
the FMFs more comprehensible through the results of 
the Ref. [35] for a Zjy rotor model as well as those in 
Ref. [33]. 
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Note added: By the completion of this work, we became 
aware that the authors of Refs. [39H4T] have indepen- 
dently achieved similar results for the fractionalized Ma- 
jorana edge states. 
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